Transmission through quantum networks 
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We propose a simple formalism to calculate the conductance of any quantum network made of 
one-dimensional quantum wires. We apply this method to analyze, for two periodic systems, the 
modulation of this conductance with respect to the magnetic field. We also study the influence of 
an elastic disorder on the periodicity of the AB oscillations, and we show that a recently proposed 
localization mechanism induced by the magnetic field resists to such a perturbation. Finally, we 
discuss the relevance of this approach for the understanding of a recent experiment on GaAs/GaAlAs 
networks. 
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It is well known that quantum transport exhibits de- 
viations from classical transport, resulting in corrections 
to the classical addition rules of conductances or resis- 
tances. A spectacular example is the Aharonov-Bohm 
(AB) effect, where the conductance of a ring is a peri- 
odic function of the magnetic flux (f> through its opening, 
with period (f>o — h/e. Since the first observation of this 
effect in condensed matter jy, many papers have been 
devoted to the study of coherence effects in transport, 
especially in the ring geometry. A first approach uses 
the Landauer formalism in which the conductance is pro- 
portional to the transmission coefficient. In this frame- 
work, disorder effects have been considered in single- 
channel ||] and multi-channel rings ||. On the other 
hand, the conductance of diffusive systems has been also 
extensively studied within the Kubo approach, where the 
weak-localization correction is related to the modulation 
by the magnetic field of the return probability of a diffu- 
sive particule [[|. Although being a transport property, 
this correction is a spectral quantity, since it is related to 
the spectrum of the diffusion equation, more precisely to 
its spectral determinant j^]. 

In this paper, we focus on the transmission properties 
of quantum networks, generalizing the original works of 
the 80's |^||. This work is motivated by recent conduc- 
tance measurements of normal metallic networks etched 
on a 2D GaAs/GaAlAs electron gas j7j. Remarkably, for 
the particular T$ network shown in fig. |l|, the magnetore- 
sistance presents large ^-periodic oscillations which are 
barely visible for a more conventional geometry like the 
square lattice. This is the first time that <f>o oscillations 
are observed in a macroscopic system where, in princi- 
ple, ensemble average due to a finite coherence length is 
expected to destroy them. 

The experimental study of the T3 network has been mo- 
tivated by the recent prediction of a new type of magnetic 
field induced localization. Indeed, it has been shown, in 
a tight-binding approach, that when the flux (f> per el- 
ementary plaquette equals </>o/2 (half- flux), the electron 
motion is completely confined inside the so-called AB 
cages H resulting from a subtle quantum interference 



effect. This surprising phenomenon has first been exper- 
imentally observed in superconducting (7^) networks [Bj, 
where it was found that the critical current almost van- 
ishes at <j) = <po/2. The standard mapping between the 
Ginzburg-Landau theory and the tight-binding problem 
Jl0| actually allows one to relate this current to the en- 
ergy band curvature, predicting a zero critical current at 
half-flux. However, it is interesting to know whether this 
localization effect still exists in normal metallic networks 
and if it could be at the origin of the oscillations discussed 
above jjj. 

The aim of this paper is threefold. Firstly, we provide a 
simple formalism allowing to calculate the transmission 
coefficient of any network made up of one-dimensional 
wires. Secondly, we concentrate on two regular struc- 
tures, the square and the T3 networks and study the flux 
dependence of the transmission coefficient which is remi- 
niscent of the butterfly-like structure of the tight-binding 
spectrum. We then consider the influence of elastic disor- 
der that we model by a distribution of the wire lengths. 
We show that the % network exhibits </>o-periodic os- 
cillations which are robust with respect to disorder and 
which are much larger than those observed in the square 
network. We also discuss the crossover from a ballistic 
(in the pure case) to a disorder dominated behaviour, 
revealed by the emergence of ^o/2-periodic oscillations 
reminiscent of the weak localization regime. This model 
gives a strong support to the interpretation of the above- 
mentionned experiment in terms of the AB cages. 

We consider a graph made up of N nodes and con- 
nected to Ni n wires (also called channels) defining the 
input reservoir and to N out wires defining the outgoing 
reservoir (see fig. [l]) . In the Landauer approach, the two- 
terminal conductance is proportional to the total trans- 
mission coefficient defined by : 

t = (!) 

where i 6 [l,iVj n ] denotes the i th input channel and 
j £ [N — N out + 1,N] denotes the j th output channel. 
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This coefficient is the sum of each individual transmis- 
sion coefficient obtained by injecting a wavepacket in the 
i th channel. We emphasize that actually eq. (|l|) assumes 
that there is no phase relationship between the different 
input channels 



FIG. 1. A piece of the T$ network (right) and of the square 
lattice (left). Black (respectively grey) dots represent the con- 
nections to the input (respectively output) channels. The 
central black dot is the input channel chosen for the bulk 
injection. 

Let us consider an incoming wavefunction in the i 
channel defined by : 

i } {x) = e- ik * + r ii e ikx , (2) 

where is the reflexion coefficient in this wire. We 
need to determine the transmission coefficient Uj giv- 
ing the probability for the wave packet to outgo into the 
j th channel. Therefore, we first solve the Schrodinger 
equation on each bond whose extremities are denoted by 
a and (3. The corresponding eigenfunctions are simply 
given by : 

i>ap{x) = -t— n — [ipa sin k(l a p ~ x ) + e^^ippsinkx] , 

(3) 

where 7 Q /3 = (27r/^>o) A.dl is the circulation of the 
vector potential between a and [3, k is the wave vector 
related to the eigenenergy by : E(k) = h 2 k 2 /2m, x is the 
distance measured from the node a and l a p is the length 
of the bond (a, /?). The current conservation on each 
internal node of the network (not connected to reservoirs) 
is satisfied if : 

M aa ij; a + M cxf}ip0 = 0> ( 4 ) 
(«,/3> 

where M is a (N x N) matrix whose elements are given 
by: 

M aa = > cot kl a /3 , M a f3 = — —- — . (5) 
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The symbol (a, (3) indicates that the sums extend to all 
the nodes (3 connected to the node a. In addition, the 
off-diagonal element M a p is non zero only if the nodes a 
and (3 are connected by a bond. Consider now the case 
where the current is injected in the channel i € [l,Ni n ]. 
The current conservation at this node writes : 

(6) 



(7) 

Finally for i 6 [I, N in ] and j e [N-N out + l, N], the con- 
tinuity of the wavefunction reads : ipi = 1 + ru and ipj = 
tij. The equations (||J|,@) constitute a (N x N) linear 
system || from which t y = tpj (j E [N - N out + 1, N]) 
can be calculated. The total transmission coefficient is 
finally obtained from eq. ([I]) by considering the AT in input 
channels. 

We now apply this formalism to the case of regular net- 
works where all the bonds have identical length I so that 
the transmission coefficient T(fc, /) is a periodic function 
of the wave vector k with period 2w/l and a periodic 
function of the reduced flux / = 4>/<po with period 1. 
In principle, the fc-dependence of the transmission coef- 
ficient can be probed experimentally if the wave vector 
k is well defined, i.e. if the energy of injected electrons 
is well controlled. Several factors like finite temperature 
or finite bias contribute to broaden this energy. This 
can be taken into account by giving a finite width Ak 
to the Fermi wave vector of the incoming wave packet. 
For example, in ref. the conductance of a single ring 
was measured and it was found that the phase of the AB 
oscillations could be varied by tuning the gate voltage, 
and thus the Fermi energy. One may therefore conclude 
that the width Afc is smaller than the period 2-k/I. These 
oscillations are very well described by a Landauer single- 
channel formalism, assuming that the ringis assymctric, 
i.e. the two arms have a different length PJl^|. 

For a given k, the flux dependence of T(k, f) has a rich 
structure which is reminiscent of the complexity of the as- 
sociated tight-binding spectrum. Here, for simplicity, we 
have chosen to average the transmission coefficient over a 
period k £ [0, 2-k/I]. The flux dependence of the average 
transmission (T(k,f))k is shown in fig. |^ for the square 
and 73 networks. One clearly observes a few peaks in the 
transmission for particular values of the reduced flux : 
/ = 1/2, 1/3 for the square lattice and / = 1/3, 1/6 for 
the 73 lattice. One can simply understand this structure 
by invoking the extended nature of the corresponding 
eigenstates that are Bloch waves with a spatial period 
proportional to the denominator of / [^3|. Due to the 
existence of the AB cages, the transmission coefficient is 




Muipi + 2J Mipipp = i(l - ru). 

For each node j £ [N — N out + 1,N], one also has 
M n^3 + M 3^P = -Mir 
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minimum at / = 1/2 for the %, network but, surprisingly, 
it is not exactly zero. This is due to the existence of dis- 
persive edges states |H| that are able to carry current 
even for / = 1/2. Therefore, T converges toward a finite 
value for the T3 network when the system size (and iVj„) 
increases, whereas T ~ Ni n in the square lattice. How- 
ever, when one injects current in the bulk of the sample, 
the transmission completely vanishes for this flux, (see 
inset fig. ||). This study shows that the cage effect, orig- 
inally predicted in a tight-binding model, also arises in a 
73 network made up of one-dimensional ballistic wires. 
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FIG. 2. Averaged transmission coefficient (T(fe, f))k/Ni n 
as a function of the reduced flux for a (8 x 8) square lattice 
(square) and a piece of the Ijj network (triangle) with 75 sites. 
Input and output channels are connected as displayed in fig. |l| 
Inset : Averaged transmission coefficient for the T3 network 
with one input channel at the center of the network. 

We now consider the case of disordered networks, the 
motivation being to see whether the cage phenomenon 
persists in such a situation. Disorder can be introduced 
in several ways (randomly distributed pointlike scatter- 
ers, or more generally, random elastic scattering matrix 
along the bonds). Here, in order to simulate random 
phase shifts on each bond, we consider a geometrical dis- 
order defined by a random modulation of the wire lengths 
while keeping the same connectivity. Denoting by Al the 
amplitude of the length fluctuations, the relevant dimen- 
sionless parameter to characterize the strength of the dis- 
order is the quantity kAl and thus explicitly depends on 
the energy. Note that the incommensurability between 
the different lengths breaks the periodicity of T with re- 
spect to k. This type of disorder also provides a distri- 
bution of areas of width 21 Al so that the oscillations are 
expected to disappear after about 1/ Al periods. In the 
following, we will focus on situations where Al/l <C 1 
and ! > 1 so that the case kAl ~ 1 may be reached 
without a sizeable dispersion of the areas. Thus, we will 
not modify the bond lengths in the phase factor e* 7 "' 3 
and the periodicity with respect to the reduced flux will 
be conserved. 



For a given realisation of disorder, T(k, f) exhibits a 
(/>o-periodic complex structure which is a signature of the 
interference pattern through the network. In particu- 
lar, the transmission extremely sensitive to k. However, 
experimentally, there is always a finite phase coherence 
length Ltp. Therefore, a two-dimensional network of typi- 
cal linear size L must be considered as a set of (L/L^) 2 re- 
gions without phase relationship. This provides a natural 
averaging mechanism over disorder realisations. Thus, 
we have chosen to study the disorder averaged transmis- 
sion coefficient {T(k, f))dis. whose variations versus the 
reduced flux are displayed in fig. ^ for fixed k and disor- 
der strength. It is clearly seen that for the square net- 
work, the periodicity of (T(k, f))dis. with respect to the 
magnetic flux is no longer </>o but 4>o/2. The (^-periodic 
oscillations have been washed out since they do not have 
a given phase. By contrast, the 0o/2-periodic oscillations 
are still present since they are related to phase coherent 
pairs of time-reversed trajectories according to the weak- 
localization picture. For the T3 network, the transmission 
coefficient remains 0o _ periodic with a large amplitude. 
This strongly suggests that the cage effect (which locks 
the phase of the oscillations) survives for this strength of 
disorder. 
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FIG. 3. Transmission coefficient {T(k, f))dis./Ni n averaged 
over 50 configurations of disorder for kl — n/3 and kAl = 1.47 
as a function of the reduced flux. 

For a finer analysis, it is interesting to compute the 
discrete Fourier transform of T defined by : 

-. 71— 1 

f(k,cj) = - VT(fc,j/n)e i2 ^> we [0,n-l], (8) 

where n is the number of sampled values of /. Figure 
D displays (\T(k,l)\)dis. as a function of the disorder 
strength kAl for different values of k. It shows that, 
when disorder is increased, (\T(k, l)\)dis. persists much 
longer for the T$ network than for the square network. 
We are thus led to conclude that the cage effect is robust 
with respect to disorder. Note that for weak disorder, 
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(\T(k, l)\)dis. depends on k but this dependence vanishes 
for k Al > 2. We strongly believe that this result ex- 
plains why a 0o-periodic conductance is observed experi- 
mentally for the "Tjj network while it is not for the square 
lattice 0. 

The behaviour of (\T(k,2)\)di S . is shown in fig. ||. It 
is interesting to see that this harmonic becomes quickly 
dominant for both networks and remains constant for 
kAl > 2. The value of this constant depends on the 
system size and converges to zero for the infinite lattice. 
Nevertheless, we leave the precise analysis of this scaling 
for further studies. 

Finally, it should be stressed that, experimentally, 
{\T(k, 2)\)dis. is further reduced by a factor e ~ 2L / L >p due 
to a finite coherence length L v , while the 4>q contribution 
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is only reduced by a factor e 
of an elementary plaquette 
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FIG. 4. Variation of {\T(k, l)\)dis./Nin versus disorder for 
different values of k after averaging over 50 configurations. 

In conclusion, we have used a simple and general 
formalism to calculate the transmission coefficient of 
any network made up of single-channel quantum wires. 
This coefficient can be simply expressed in terms of a 
connectivity-like matrix. We have used this formalism to 
study the AB cage phenomenon in the T3 network and 
we have shown that this effect is robust to a moderate 
amount of elastic disorder. As a consequence, the AB 
oscillations with period <po persist in the infinite T3 net- 
works whereas they vanish in the square lattice. 

We acknowledge H. Bouchiat, P. Butaud, R. Deblock, 
G. Faini, D. Mailly, R. Mosseri, C. Naud and B. Reulet 
for useful discussions. 
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FIG. 5. Amplitude of the second har- 
monic {\T(k,2)\)dis. /Ni n versus disorder for different values 
of k after averaging over 50 configurations. 
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